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_ , ' Abstract: 

^ ' The dynamics of the driven single band tight binding model for Wannier-Stark systems is for- 

. mulated and solved using a dynamical algebra. This Lie algebraic approach is very convenient 

I for evaluating matrix elements and expectation values. A classicalization of the tight binding 

' model is discussed as well as some illustrating examples of Bloch oscillations and dynamical 

■ localization effects. It is also shown that a dynamical invariant can be constructed. 

o 
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1. Introduction 

The celebrated single band tight binding system 



k> , ^ +00 +00 

^; ^ ^ ~^ X] + 1| + 1"- + X] (1) 



models a space periodic system with period d in a (possibly time dependent) linear field. Here, 
n numbers the sites and |n) are the Wannier states with (n|n') = 5nn' ■ In © only nearest 
neighbor interactions are taken into account. In this model, the periodic field free system has 
only a single band with dispersion relation 

E{k) = cos (Kd) (2) 

where k is the Bloch index and A is the band width. In the simplest case, the field F is constant, 
a dc-field. More complicated is the combined ac-dc-system with time periodic driving; an often 
considered case is the harmonic driving, F{t) = Fq — Fi cos(u;t) . 

The dynamics of the driven tight binding system is quite involved and, despite of the large 
number of previous studies, of increasing interest, in particular in view of the recent progress in 
studies of the dynamics of ultracold atoms in standing wave laser fields. For more information, 
see |ll2j and the references given there. It is well known that the tight binding system allows 
an analytic treatment and various approaches have been proposed (for early studies see |.'-{|4|5j ) . 
Quite generally, however, the derivations are quite tedious. Here, we recommend a treatment 
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based on the dynamical Lie algebra |6l7l8j which appears to be favorable because of its generality 
and simplicity. This approach allows a straightforward evaluation of the time evolution operator, 
matrix elements, expectation values and dynamical invariants by purely algebraic operations. 

2. The algebra 

The three operators N, K, where N is hermitian and K unitary with commutation 
relations 

[k,N]=k, [klN]=-k^, [k^,k]=o (3) 

form a closed Lie algebra C. This shift-operator algebra [9 is obviously different from the ubi- 
quitous oscillator algebra {n, d, }, but some features are similar (see also |lflj for a discussion 
of the more general quantum boson algebra which contains both algebras as limiting cases). 
The dynamics generated by the hermitian Hamiltonian 

H = G{t) {k + k^) + F{t)N (4) 

with real valued, possibly time dependent functions F and G can be conveniently studied by 
algebraic techniques. 

A realization is the tight binding model with 

+00 +00 +00 

iV = ^ n\n){n\ , ^ = \n){n + l\. , K"^ = ^ |n + l)(n|. (5) 

n=—co n=~oo n=~oo 

It should be noted, however, that this algebra appears also in different context Pj and therefore 
some general considerations seem to be appropriate. First, one can easily show that k and k^ 
act on the eigenstates of N, N\n) = n\n), as shift- or ladder-operators: 

k\n) = \n-l), k^^ln) = \n+l) , (6) 

which fixes the eigenvalues of at no + n with n E Z up to an arbitrary value of no in the unit 
interval. This value can be fixed if one considers the algebra as a subalgebra of a bigger one 
by adding an antiunitary operator representing time inversion. This leads to the two possible 
cases of no = (bosonic) or no = 1/2 (fermionic). See [TTj and references given there for more 
details. Here we are interested in the bosonic case, i.e. 

N\n) = n\n) , n G Z. (7) 

In context of the tight binding system, is a 'position operator': the expectation value (A^) = 
(TplNltp) is the mean position on the lattice and pn = \{n\'ip)\'^ is the population probability 
of the 'lattice site' at position n. This is the physical system we have in mind. It should be 
noted, however, that the same algebra appears in different contexts as, e.g., for the plane rotor 
[IJ±]=±J±,[J+,J_] = 0. 
The eigenvectors of K with eigenvalues e"^ are \k) = J2n=-oo e^"'^!"-) '■ 

+00 +00 
k\K) = e'"'^|n - 1) = ^ e'^ X] e''"''|m) = e''' \k) . (8) 

71=— 00 m=— 00 

These 'Bloch states' are 2tt periodic and normalized as 

+00 

{k\k') = X S{k - k' - 27m) = 52n{K - k') (9) 
n= — 00 
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where b2-K is the 27r-periodic comb function. 

The representation of the operator N in this basis is 

(K|iV|K') =(52^(K-K')i^- (10) 

The algebra C = {K, ,N} has the radical TZ = {K, K^} , the simple part <S = {N}, and can 
be decomposed into the semidirect sum C = IZ-^S as, e.g., described in For a Hamiltonian 
H = Hg + Hji with G IZ and G 5, the time evolution operator can be factorized: 

U = UsUr (11) 

with 

ih^ = HsUs, ih^={Us'HRUs)UR. (12) 
(Moreover, Us can be factorized into simple parts if S is only semisimple; for more details see 

IH1-) 

This product decomposition has several advantages in comparison with the pure exponential 
solution, in particular it provides a global solution [12ll3j and the calculation of expectation 
values and matrix elements is simplified as will become clear later on. 

The initial step of any application is the evaluation of all necessary T-evolved operators A of 
interest, i.e. 

g^adf^ := e'^Ae-'^ 

2 3 

= i + z[f,i]+|j [f,[f,i]]+|j [f, [f,[f,i]]]+ ... (13) 

with z G C for all F G Trivially we have e^'^'^^F = F . Here we certainly need the evolved 
operators of our algebra C = {K , , N} . Because K and commute, we have 

e^^^^F(i^t) = F{K^) , e^''^^V(J^) = F{K) (14) 

and the nontrivial expressions are 

-,zadKTCT_ tCt I „r> „zadiVr> -z t> „zaxlK'l tCt _ tCt _ r>t zaxiN 



e 



N = N + zk , e'^'^^k = e-'K , e'"""^^ N = N - zK^ , e'^"^^ = (15) 



which can be easily obtained from (|13() using ©. 

3. Time evolution operator 

For the tight binding Hamiltonian Q, where for simplicity we introduce the notation gt = 
G{t)/h, ft = F{t)/h, the simple part of the time evolution is lUs = ft^^Us with solution 

Us{t) = e-"^^^ , r^t= /Vrdr (16) 

Jo 

and the remaining equation of motion i Un = {Ug^Hn Us) Ur for the radical part can be solved 
in a second step. Using the relation (jlSf) we find 

Us'Hr Us = e'^^'^'^^HR = hgt {e^'^^'k + e+'^^^k^) (17) 
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and therefore 

Udt) = e-'(>^*^+>^*^') , xt = fgr e-''-'^ dr (18) 

Jo 

and finally 

U{t) = Us{t) UR{t) = e-i^'^e-'>^'^e-'>^*-^^ , (19) 

the Wei-Norman product form of the time evolution operator j8ll2ll3j which is, in fact, a version 
of the so-called momentum gauge in this case. 

Furthermore, a series expansion in powers of the ladder operator will be useful, which can be 
obtained by means of the generating function for the Bessel functions 

+0O 

J„(2u)B". (20) 

n=— oo 

Identifying B = e~^^'^*^'^^'^^ K where xt = Ixtl^"^*^* , equation H18|) can be rewritten as 

+00 

UR{t)= Jn{2\xt\)e-'<'f'^+-/^^K^. (21) 



In various applications, matrix elements of the time evolution operator are required. Making 
use of 

(K|e-™*|K') = y (K|n)e-'""(n|K') = e'^^'^'"'^-") = 62^1^' - k - u) , (22) 

n n 

the matrix elements in the Bloch wave basis can be directly read off from (|19)): 

{K,\ij{t)W) = 62^{n' - ,^-r^^)e-'^^xt\cos(K'-.j>t) _ (23) 

Matrix elements of the propagator (|19() in the basis |n) follow immediately from H21() and the 
ladder property K'^\n') = \n' — n): 

Unn'{t) = e-i(-'-)(^.+-/2)-in,. j^,_^(2 \xt\) (24) 

which coincides, of course, with the result derived many years ago by Dunlap and Kenkre |^). 

For completeness, we should also state the explicit results for the most frequently studied 
cases: 

(1) For time independent functions gt = go and ft = fo the integrals in (fT^ and ((TH|) yield 

Vt = fot, Xt = ^ e-i^o*/2 sin(/ot/2) . (25) 
Jo 

Note that at time t = Tb = 27r//o the evolution operator is equal to the identity 

IJ{Tb) = e-^2"^ = / (26) 

because of e"'^'^^ \n) = e"'^'^" |n) = \n) . Therefore the dynamics is periodic with the Bloch 
period Tb and Bloch frequency ujb = fo- It is also of interest to compare the product form (|19|) 
of the propagator with the pure exponential one which is trivial in this case, namely 

IJ(^^^ ^ Q-iHt/h ^ g-i {goK+goK''+foN)t _ (27) 
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The non-obvious identity between ()27|) and ()19l) becomes clear in view of the generahzed Baker- 
Campbell-Hausdorff formula P| 

^c{X+f}Y) ^ g/3(e--l)y ^aX ^ ^aX g/3(l-e-")y ^28) 

for shift-operators X, Y with commutator [X, y] = y . 

(2) For harmonic driving, 

ft = fo- fi cos(cjt) , gt = go, (29) 

we have 

Vt = fot sin(wt) (30) 

and - using again the Bessel expansion (^0]) - 

ft ^ ft 

Xt = go dre-i^o^+^^^^-(-^)=go Yl •^-(^) / dre-^^"^ 

oo ^ 

= 250 5Z '^-i^)— e"' sin(c^,t/2) , c^, = - i^o; 7^ . (31) 

1'=— 00 

This is an oscillating function of time. For resonant driving, 

ujb = nu) , n = 1, 2, . . . , (32) 

the integration of the nth term in the sum (|31j) yields a linearly growing term, which dominates 
the oscillating rest of the sum for long times, i.e. we have 

Xt-lnt/2, Jn = 2goJn{^)- (33) 
Later on, in section |H1 some consequences of this resonant behaviour will be discussed. 

(3) The general case of a combined dc- ac-system can be treated in a similar manner. Let us 
consider the case 

ft = fo + ft , ft+T = ft , gt+T = gt (34) 

with /o chosen according to ftdt = . Again we consider the case of resonant driving, 
T = uTb, with Tb = 2tt /ujb , = fo ■ Fourier expansion of the periodic part of the force 

+00 

ft= ' ^0 = 0' (35) 

/x=— 00 

with UJ = 2tt/T yields 

m = fot + Y^[ e"^' - 1] = iOBt + Vt (36) 
where fjt is T-per iodic. A second Fourier expansion 

-00 



gte-'^' 



V a, e''^-* , a. = 2^ / gt e-'^-*-^* dt , (37) 



allows the evaluation of xt '■ 

rt rt 



Xt= [ gre-'^^dT= [ 5,e-^-^^-^'^^dr = a„t + V:^[l-e---*] (38) 
Jo Jo I'^v 
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with 

uji, = ujb — ■ (39) 
Note that this is again a sum of a hnear growing and a T-periodic part: 

Xt = ant + Xt , Xt+T = Xt ■ (40) 

The coefficients of two Fourier expansions (|35j) and (|ci7|l are, of course, not unrelated for constant 
go- Let us confine ourselves here for simplicity to the case of a symmetric resonant driving: 

oo 

/t = /o + /m cos(mwt) , UJB = nLO, (41) 

m=l 

with - inserting fo = - 

oo 

r]t = uJBt + ^ /3m sin(ma;t) , Prn = fm/muj . (42) 

m=l 

Using now the generating function for the infinite-variable Bessel functions |14j (for a recent 
application of these little known functions to a two-dimensional tight binding system see ^H] ) , 

oo +0O 

exp (i l^ra Sin mn) = ■Wm}) e'"" , (43) 

m=l u=—oo 

the remaining integral for xt can be evaluated analytically: 

ft 



Xt = go 



f e-'^-dT = goJn{{(3rn})t + 2goy^M{(3^})—e-''^^'/^sm{uj,t/2). (44) 



Similar to single frequency driving, dynamical localization effects (see section|SI) may be observed 
for system parameters leading to a zero of the infinite-order Bessel function Jni{Pm})- This 
deserves future studies. 

4. Quasienergies 

In the case of a combined dc- and time periodic ac-driving 1)34(1 under resonance conditions 
T = uTb, the time evolution operator (|19|) over a period T simplifies. Inserting (|36|) and (|4U() 
with fjx = fio = and XT = Xo = 0) we find 

U{T) = q-'^VtN ^-\xtK ^-ix'tK^ = ^-lar^Tk ^-ialTk^ ^45) 

which commutes with K and allows the construction of simultaneous eigenstates 

i^|^,(r)) = e-|^,(r)). (46) 

IJ{T) |V^«(T)) = e-ia„Te+'^-ia;re-i'= ^ ^-i..T _ (47) 

The quasienergies are identified as 

= ane+''' + a^e-^" = 2|a„| cos(«: + 93) , (48) 
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the dispersion relation for the quasienergy. Here 4:\an\ is the width of the quasienergy band and 
(p, the phase of the Fourier coefficient an = \an \ e^'^, is zero if gt and ft are symmetric in time. 

It is also of interest to construct explicitly the time dependent quasienergy (or Floquet) states. 
As can be easily seen, the states 

\Mt)) = U{t) 1^) = -1= ^ei«'''-ife'^+"'+x^-'^)|n) , Kt = K-rjt (49) 

n 

are solutions of the time dependent Schrodinger equation 

d 1 

{i—--H)\i;,{t)) = (50) 

and, simultaneously, eigenstates of K : 

k\ij^it))=e'^^\Mt)), (51) 
the so-called Houston states. Using again and (flUl) '). we have 

n 

Therefore the state 

Kit)) = e+'^«* \Mt)) (53) 
with given in 1)48^ is T-periodic, \Uf^{t + T)) = \ui^{t)) and solves 

{i^-^H)\u^(t))=eKit)), (54) 

i.e. it is a quasienergy state and e is the quasienergy. From (|49|) we see that the Floquet states 
extend over the whole lattice. As a final remark, we also note the obvious identity \uk{T)) = \k) . 

5. Expectation values 

The time dependence of expectation values follows immediately ||8j from the relations H14|) - 

and therefore 

{k)t = e-'''* (K)o = e'('^-^*) \K\ (56) 

with (i')o = K = \K\ e'*^ . From 

k^{t) = e-2i^* k^ and {k'^)t = e'^i^* {k^)^ , (57) 
we see that, up to a phase factor, the variance is constant: 

^l{t) = Uk%-{k)l\ = ^l{Q). (58) 
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The time dependence of the position operator is a bit more interesting: 

N{t) = e'^'^^ e'^'^ N e-'^'^ e''^'^^ = e'^'^^ {N + ixtK) e-'^*^^^ 

= N + i{xtK-x*tK^) (59) 
and therefore, using xt = Ixtls"^*^', 

{N)t = {N)o + i{xt{K)o-X*t{K%) = {N)o + 2\K\\xt\M(Pt-i^). (60) 
Introducing the anti-commutator J = NK + KN = [N , K]j^ , the time evolution of iV^ is 

N\t) = [N + \ ixtK - x*tK^)f = iVVi ( XtJ-xV^)-xlK^-xfK^^ +2 \xt? (61) 
and, with 

(J)o = J= |J|e'^ , {k% = L = \LW\ (62) 
the expectation value evolves as 

{N^)t = {N^)o + 2\J\\xt\ smicPt - ^) + 2 |xt|'(l - \L\ cos(20t - i^)) . (63) 
Finally, the time evolution of the variance is given by 

A^(t) = {N%-{N)^ 

= A%{0) + 2\xtmi-\L\cos{2(t>t-i^)-2\K\Hm'^{^t-K)} 

+2 \xt\{ 2{N)o\K\ sin((/>t - k) + \J\ sin(</.t - ^) } . (64) 

The dynamics of the expectation values therefore depends on three complex coherence parame- 
ters which are explicitly 

= ' 'J = ^C^^ - (^n-lCn , L = ^c;_2C„ (65) 

n n n 

if the initial normalized state is specified as jV') = Yln^ri \'n). 

Sometimes it may be more convenient to replace the unitary shift operators K and K'^ by the 
hermitian operators C and 5, 

k = C + iS, k^ = C-iS, (66) 

whose expectation values allow a direct interpretation. Clearly, also these operators commute 
and the commutators with the position operator are 

[C,N]=iS, [S,N]=-[C. (67) 

Rewriting xt defined in equation (|18() as 2xt = ut — ivt with 

ut = 2 Qr COS r/r dr , vt = 2 Qr sin r]T dr (68) 
Jo Jo 

the solution for the position operator in and (|^n|) is 

N{t) =N + vtC - utS (69) 
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{N)t = {N)o + vt{C)o - ut{S)o (70) 

and, with J = [N, C]-^- + i [N, S]-^- , the time evolution of iV^ in and is rewritten as 

Ar2(i) = N"^ + vt[N,C]+-ut[N,S]+ + v^C'^ + u'is'^ -2utVtCS (71) 

(iV2)i = {N\ + [iV, C']+)o - ut{ [N, S]+ )o + v^{C\ + 7x?(52)o - 2utVt{CS)o (72) 

with {C'^)o + (5'^)o = 1 • The time evolution of the variance of the position can then be 
formulated in the convenient form 

A%{t) = A%{0) + 2vt - 'int + v''^ A^c + - 2utvt A^s (73) 

where is the covariance of the expectation values of the operators A and B at time t = 0: 

Ai^ = (l[i,B]+)o-(i)o(^)o. (74) 

Note that the relations for the expectation values derived above are valid for pure states as well 
as for mixed states (see |16|17|18fl9] for an application of the tight binding system using density 
matrices) . 

6. Classicalization 

Recently, a classicalization of the tight binding model with Hamiltonian 

H = 2G{t) cos{p6) + F{t)q/d . (75) 

has been discussed |2l22j : related observations can also be found in jl8l23j . In this classicaliza- 
tion, the operators N and C = {K+K^)/2 are replaced by phase space functions. The parameter 
6 = d/h depends explicitly on h which implies, of course, that this 'classicalization' differs from 
the usual classical limit of quantum dynamics. It has been observed that the classical dynamics 
generated by ((751) . 

r)J-f f)J-f 

'* = -W = -d- ^ = ^ = -^assHpS). (70) 

show a surprising agreement with the quantum one. Here we will analyze this correspondence 
from an algebraic point of view. 

First, we can again generalize and consider a classical Hamiltonian (|75|) with time dependent 
coefficients. Introducing the dimensionless phase space functions C{p) = cos{p5), S{p) = sm{p5) 
and N{q) = q/d, with Poisson brackets 

{C,N] = +\S{p) , {S,N} = -lC[p) , {C,5} = 0, (77) 

the set {N, C, S} forms a closed Lie algebra with Lie bracket { , } and therefore the dynamics 
induced by the Hamiltonian H{p,q,t) = G{t) C{p) + F{t) N{q) can be evaluated again by purely 
algebraic techniques. 

Moreover the classical algebra and the quantum algebra studied above are isomorphic in the 
present case which is evident from the mapping 

N < — >N , C i — >C , Si — >S (78) 

where the operators C and S are defined in H66() (the Lie brackets map according to { , } < — > 
[ , ] ) . This implies the equality of the dynamical evolution of quantum operators and clas- 
sical phase space functions in this case. In particular, the evolution of the expectation values 



10 



agrees. There are, however, some differences, e.g. for the mitial conditions. Whereas there is no 
hmitation in the classical case, the quantum covariances are limited by uncertainty relations. 
Furthermore, the equations of motion 

C = fS, S = -fC, N = -gS, (79) 

with / = F/h, g = G/h, are linear which implies that the classical dynamics is regular, i.e. not 
chaotic. 



7. The dynamical invariant 

The driven tight binding system Q) possesses a dynamical invariant, very similar to the 
harmonic oscillator with time dependent frequency, where the so-called Lewis invariant / |2UI21j , 
a time dependent constant of motion, 

plays an important role. 

The dynamical algebra offers a convenient technique for constructing such an invariant |6|7j . 
Writing the invariant as a linear combination of the basis of the algebra, / = AjFj , as well 
as the Hamiltonian and inserting these expressions into equation (|8fl|) using the commutator 
relations, one obtains a set of linear first order differential equations for the coefficients Xj{t). 
In the present case, we write 

i = jN + xk + x*k^ (81) 

with 7 G M because the invariant can be chosen to be hermitian. The procedure described above 
|6l7j leads directly to the differential equations 

X = iiftX-ga), 7 = 0, (82) 

i.e. 7 is an arbitrary scaling constant which can be chosen as 7 = 1 and the solution for 7^ is 

Xt = -ie'^'Xt (83) 

were ijt and Xt are defined in equations (|16j) and (|18j) . The full expression for the invariant is 
therefore 

i{t) = N{t) + Xtk{t) + x;k^t) = /(o) = n (84) 

where the time evolved operators K{t) and N{t) have already been calculated in equations H55() 
and (|59j) . The time evolution conserves the commutator relations © and therefore K{t) and 
K^{t) are still acting as ladder operators on the time dependent eigenstates \n,t) of N(t) with 
the (time independent) eigenvalues n S Z. Following Lewis and Riesenfeld |U, one can also 
construct the general time evolution from the eigenstates of the invariant. 

In view of the classical version of the tight binding model in the preceding section, we will 
also express the dynamical invariant in terms of the hermitian operators C and S given in 
which yields after some algebra 

I{t) = N{t) + {ut sin rjt - vt cos r]t) C{t) + {ut cos rjt + vt sin rjt) S (t) . (85) 
In the classical version (jTSJ, the invariant is the phase space function 

I{p, ^' ^) ^ ^ ®™ ^* ~ ''-'^ cos(p5) + {ut cos rjt + Vt sin r]t) sm{p6) . (86) 
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where p = pt and q = qt evolve under the Hamiltonian eqautions of motion H76|) and 

Iip,q,t)=I{po,qo,0) = ^ (87) 

is a constant of motion. 

8. Oscillating versus breathing modes and dynamic localization 

Despite of the algebraic simplicity, the tight binding dynamics shows some non-intuitive fea- 
tures, even in the case where the Hamiltonian does not explicitly depend on time, where the 
famous Bloch oscillations are observed. Much more phenomena can be found for driven system, 
as for instance dynamical localization effects j4l23llj . A discussion of these phenomena is far 
beyond the scope of the present article. Some quite general features, however, can be directly 
seen from the dynamics of the expectation values for the position operator and its variance. 
For simplicity we will use the classical description of the tight binding dynamics outlined in the 
preceding section. To avoid misinterpretations, it should be recalled that the Bloch oscillation 
remains, of course, a pure quantum phenomenon because the 'classical' system is /i-dependent. 

First, the general dynamical behaviour is strongly influenced by the initial distribution, more 
precisely by the expectation values of C{p) = cos{p5), C^ip) = cos^(p(5), . . . which satisfy the 
obvious bounds -1 < (C), {S) < +1 and < (C^), (5^) < +1 . We will assume in the following 
that the initial classical phase space distribution is symmetric in the position g, which implies 
(iV)o = and A^^ = = 0, i.e. equations ((ZQl) and read 

{N)t = vt{C)o-ut{S)o (88) 

(t) = A%iO)+ Ale + Als - ^utvt Als . (89) 

Now the dynamics is most strongly influenced by the localization properties in the momentum. 
Let us distinguish two extreme cases: 

(1) If the initial distribution is sharply localized in the vicinity of momentum po, we have 
(C)o ~ cos(po<^)) ('S')o ~ sm{pQ5) and A^^-, A|.^ ^ ^cs ~ ^ therefore 

{N)t ^ vt cos{po5) - ut sin(po5) , (90) 

A^(t)«A^(0), (91) 
and the distribution moves in space with constant width. This is an oscillatory mode. 

(2) If the momentum distribution is broad and approximately constant over a period of cos{p5), 
we have (C) « (S") and A^f^ f» A|^ r:^ 1/2 and A^^ and therefore 

(iV)i«0, (92) 

Al{t)^Al{Q) + \{vl + ul) (93) 

and the distribution is frozen in space with a time dependent width. This is a breathing mode. 
(Alternatively, equation (|93() can be derived directly from ()24() for an initial distribution localized 
on site ra [1].) 
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In addition to these influences of the initial conditions, there are, of course, various effects 
from the time dependence. For constant fields {g = go, f = /o) vt = 2gofQ^ sm{u)Bt) and 
ut = 25fo/o~^(l ~ cos{iVBt) ) are periodic in time oscillating with the Bloch period Tb- Therefore 
also {N)t oscillates with period Tb in an interval of width (-/V)max — (-^)min = 4|5o/(7^l { {^)o + 

Much more complicated is the case of an explicitly time dependent driving. Here, we will 
only briefly mention the important dynamical localization observed for resonant driving ft = 
fo — fi cos(a't) with ujb = nu , n = 1, 2, . . . , (see also the discussion following equation (PU)) ). 
Here, xt grows linearly with time (|33|) with superimposed oscillations and therefore ut ~ 7nt, 

~ 0. Then the dominant terms in equations ()88|) and H89|) are 

(iV)t = -7„(S)ot , A%{t) = A%{0)+jlAlst^. (94) 

(valid for a space symmetric initial distribution). Therefore, for large times, the width of the 
distribution increases linearly in time. This strong dispersion can, however, be suppressed by 
adjusting the field parameters to a zero of the Bessel function, i.e. 

7n = 2goJn{^)=0, (95) 

an effect known as dynamic localization |4| (see also jll23j for more details). 

9. Single band model 

The single band model is an extension of the tight binding model by replacing the cosine 
dispersion relation (j2j by a more realistic periodic function E{k). Instead of the Hamiltonian 
((IJ) one considers the generalization 

oo 

H = hJ2 {gmit) + gUt) K^"^) + F{t)N = Hr + F{t)N . (96) 

m=0 

Here, the algebra is extended to the set C = {K"^,K^"^,m G N, iV} with radical TZ = 
{K"^,K^"^,'m € N} . The subsequent analysis follows exactly the same lines as in the tight 
binding model. By means of the auxiliary relations [K"" , N] = 2"'-^k"' and [K^"" , N] = 
_2>n-ij^]m Qj^g obtains 

oo 

U-s'HrUs=Y. {9rn{t)e--^'-'^^^k-^ + gUt)e-^-^'-'^^^k^-) (97) 

m=0 

and the time evolution of the radical part of the algebra is given by 

Udt) = exp (-iE™=o {Xn.it) k"^ + X*m{t) i^t™)) (98) 

with 

Xm{t)= rdT5™(T)e-2-^''^. (99) 
Jo 

The full time evolution operator is again U(t) = Usit) UR{t) , where Usit) is still given by ((1^ . 

For most applications, however, the coefficients gm will be independent of time. In such a 
case, the dispersion relation is given by the Fourier series 

oo 

E{k) = n'£{9m e^™'^ + g*rn e"''"'' ) . (100) 

m=0 
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Matrix elements of U{t) in the Bloch basis are similar to (|23() : 

{K\U{t)\K') = -k' - r,i)e-'(^(«'*)+^*('^'*)) , (101) 

with /{k, t) = Ylm Xm{t)e^"^'^ . The matrix elements in the \n) basis are, however, more compli- 
cated than the tight binding expression ()24() . 

The time evolution of the ladder operator K{t) is still given by and, using 

e'^^K"^ N = N + z2"'-^k'^, (102) 

the evolution of the position operator is 

00 

N{t) = iV + i ^ 2™-i ( xm{t) - X*m{t) K^"") (103) 

m=0 

with expectation value 

00 

{N)t = {N)o + i J2 2"^"'(x™(t) - X*m{t) {k^n ) • (104) 

m=0 

Finally it should be noted that also in this case a classicalization is possible with classical 
Hamiltonian 

H{p,q,t)=E{p6)+Fq (105) 
as discussed above for the tight binding model. 

10. Concluding remarks 

The driven tight binding system is in many aspects very similar to the driven harmonic 
oscillator. It can be treated algebraically by means of ladder operators, there exists a dynamical 
invariant, and one observes a close correspondence between quantum and classical time evolution. 
Some of these features of the tight binding dynamics have been discussed in the present paper. 
There are, of course, still a number of interesting questions to be answered as for instance the 
role of the coherent states jl 1124125] and the relation between the invariant and the quasienergies 
for a periodically driven system |26|27j . Moreover, an extension of the algebraic technique to 
treat Bloch-Zener oscillations in doubly periodic structures jUj or the recently investigated two 
dimensional case |15|22|28p29j . Work in these directions is in progress. 
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